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DISTRIBUTED MODEL PREDICTIVE CONTROL FOR ENERGY SYSTEMS 

IN MICROGRIDS 
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Abstract. This paper presents a flexible and modular control scheme based on distributed 
model predictive control (DMPC) to achieve optimal operation of decentralized energy systems 
in smart grids. The proposed approach is used to coordinate multiple distributed energy 
resources (DERs) in a low voltage (LV) microgrid and therefore, allow virtual power plant 
(VPP) operation. A sequential and iterative DMPC formulation is shown which incorporates 
global grid targets along with the local comfort requirements and performance indices. The 
preliminary results generated by the simulation of a studied case proves the benefits of applying 
such a control scheme to a benchmark low voltage microgrid. 


1. Introduction 

The progressive shift towards decentralized generation in power distribution networks has 
rendered the problem of optimal operation of distributed energy resources (DERs) to be increas¬ 
ingly constraining. Indeed, the integration of flexible deterministic energy systems combined 
with the strong penetration of uncontrollable and stochastic renewable energy sources (RESs) 
has pushed the gird to its operating limits [T]. Micro grids represent a promising concept to 
face the latter issue; by interconnecting and thus monitoring the different DERs, the aggregated 
operation strategy enables the provision of voltage support and other ancillary services to the 
local distribution system operator (DSO) [31|31I11[S]. Nevertheless, the growing integration of 
polygeneration systems (e.g. combined heat and power units) and heat pumps raises the need of 
considering both electrical and thermal power requirements while establishing optimal operation 
strategies. 

Several studies have proposed robust control mechanisms to achieve stable operation of low- 
voltage (LV) microgrids to cope with the abovementioned objectives. From the DSO perspective, 
centralized hierarchical control schemes are commonly selected (e.g. Hi) , establishing set 
points for the DER load profiles to reach the grid operational targets. The results presented in [7] 
demonstrated the advantage of implementing the following control strategy; the voltage profiles 
at the grid buses are highly improved through optimal DERs scheduling. However, regarding the 
interest of DER owners connected to the microgrid, local strategies may highly differ from the 
DSO operating objective [6]. From the point-of-view of the end-user (e.g. a residential building), 
model predictive control (MFC) architectures are usually applied to minimize operating costs 
while satisfying the different comfort and service requirements of the smart grid actor (e.g. BM 
unmi]). Nevertheless, a lack of communication between the different MFC units while defining 
the local strategies might cause critical grid operating conditions [5]. In order to face these 
issues, authors in [12] presented a distributed model predictive control (DMFC) architecture. 
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coordinated through an independent system operator (ISO) in order to steer towards a global 
grid objective of peak-shaving. However, a robust control regarding the violation of power flow 
limit in the LV network slack bus and variable bound setting are not discussed in that study. 

This paper presents a flexible, modular and robust control architecture based on an DMPC 
problem formulation m- The fully connected, iterative and independent MFC algorithm estab¬ 
lishes the optimal set points of the controllable DERs connected to the LV microgrid. Besides 
optimizing the local control objectives, the coordination of the regulators allows the ISO to act 
as a virtual power plant (VPP) while providing ancillary services to the DSO. Indeed, by sup¬ 
plying day-ahead load predictions - computed considering a cost performance index - to the 
DSO, it is able to act with more certainty in the electricity market. The structure of the paper 
is the following: Section II explains the control scheme and Section III defines the case study 
considered. Section IV shows the advantages of using the proposed DMPC through simulation 
results. Section V finally provides concluding comments about of the proposed control method. 

2. Control Algorithm 

In addition to avoiding a conflict of interests among the different DER owners and the local 
DSO [B], the major advantage of the decentralized control remains in the computational effort 
required to solve the problem formulation [121 [E]- In fact, by splitting the large, centralized 
optimization process into multiple, smaller sub-problems, solving time and flexibility of the MPC 
is strongly improved. Nevertheless, in order to ensure that local actions are not threatening the 
global system stability, coordination between the different regulators is still required to satisfy 
specific coupling constraints. As presented in several control schemes have been proposed 
to cope with the coordination challenges. Nevertheless, no generic design approach has yet been 
developed, leading to only case-specific, tailored DMPC problem formulations. 

The DMPC architecture proposed in this paper relies on the combination of a sequential 
and iterative solving approach. As discussed in regarding the information exchanged, 
the local regulators require comprehensive models of the different subsystems. Nevertheless, 
given the global system objective considered in this study, solely the future predicted control 
variables need to be transmitted to the different controllers, without any knowledge of the other 
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Figure 1. DMPC computation scheme. 
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DER models. Since the defined coupling constraint solely relates to the net power flow at each 
bus, efforts concerning the exchange of information and controller synchronization are strongly 
improved. The following paragraph details the specific DMPC algorithm developed for the 
optimal operation of DER nnits located in an LV microgrid. 


2.1. DMPC architecture. In the presented control algorithm, the local regulators are mon¬ 
itored through an ISO which stores the predicted load profiles of the MFCs connected to the 
microgrid. The latter independent grid agent solely consists of a simple data managing system, 
representing the aggregated interests of the local network actors. To solve the DMPC problem, 
a sequential updating approach [12] is applied during which each regulator solely interacts with 
the ISO. As shown in Figure [TJ the MPC hence calls the central ISO to gather the previously 
predicted operating plans of the remaining DERs (1), computes the local strategy (2) before 
submitting it to the ISO (3). 

This process is performed in an iterative manner until the aggregated load plan variance 
perceived by the ISO reaches the defined convergence criteria e (Algorithm [1]) . In order to 
initialize the solving scheme for the first iteration (i.e. I = 1), the operation plan defined 
during the previous time step, — 1), is considered for units located later in the solving 

scheme. Each controller solves a standard MPC problem formulation. Also as shown in [TC] . 
the performance index is composed of local objectives and an aggregated input target (I). The 
thermal and electrical loads related to the local DER and service requirements are computed 
considering a I hour time step (sampling time) with a 24 hours prediction horizon. The control 
inputs Mj are defined as simply being the actual state of operation/charge of the different DERs 
within the subsystem (i.e. building) i. In order to prioritize thermal comfort requirements, a 
strong penalty cost is applied when violating the soft constraint (3). 


3. Case Study 

This case study presents the optimal operation of multiple DERs connected to an LV microgrid 
through flexible DMPC. The global objective of the distributed control architecture is to provide 
the day-ahead, aggregated load profile to the respective DSO of the LV network through the 
grid ISO. The forecasted consumption curve highly improves the DSO bargaining power on the 
daily electricity market since the future load profile is assumed to be determined a priori. In 
exchange for the transmitted information, the microgrid end-users might benefit from a reduced 
electricity tariff or similar economic incentives from the local DSO. 

To reach the microgrid forecasted load profile, the DMPC regulators are however constrained 
to respect the aggregated consumption predictions and thus, need to steer towards the defined 
profile. In order to account for the stochastic related to the specific end-user behavior, the 
aggregated consumption is allowed to vary within a predefined band in which no penalty is 
perceived by the microgrid. Nevertheless, when exceeding the upper or lower profile limit, a 
non-compliance fee (Table H]) is charged to the ISO since the economic advantage of the DSO 
has vanished. Hence, in the following case study, each grid bus (i.e. end-user) is equipped 
with an MPC regulator which optimizes the hybrid performance index (Eq. I) which includes 
the local operating expenses (Eq. 2) and comfort penalties (Eq. 3) in addition to the global 
non-compliance cost (Eq. 4): 
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{j{Xi{k),Ui{k),ai{k)) = fopex{Ut{k)) + fconf{Xi{k)) + fgrid{U^{k),a^{k)) 

(1) < Xi{k) = {xi{k + l|fc), ...,Xi{k + iV|fc)} 

\^Ui{k) = {ui{k\k), ...,Ui{k + N - l|fc)} 

(2) fopex {Udk))= Cei (k) {P+ (k)-P-ik))+ Cfuel (A:)rh+„,, (fc) 

( 3 ) fconf (Xi (fc)) = Cp^conf \fi{k)-fi{k)\ 


(4) fg„d{UPk)) = Cp,g„d PHk) - U+(fc) + E - 1) 

\ / 

where the symbol " indicates the predicted values of the corresponding state, x, and input, u, 
while +/— superscripts indicate specific input and output flows (e.g. power, P) respectively. 

Algorithm 1: DMPC algorithm 
Data: 

- R Controller set, R G N* 

- I Main loop iterator 

- U\k) Aggregated load profile of all controllers at time step k 

- cr-(fc) Aggregated load profile of controllers j, j G R — {i} 

- ul{k) Predicted load profile of controller i at time step k, i G R 
Result: Ui{k), i G R 

Initialize: I = 0, U'^{k) = 0; 
while U^{k) — U^~^{k) < e do 
^ Z + 1; 

for i = 1 to n do 
it I = 1 then 

1 — 1 n 

= E P](k) + E - 1); 

j=i i=i+l 

else 

i — 1 n 

a\{k)=Y.u^Ak)+ E u\-\k)-, 

0 = 1 j=i+l 

end 

Compute MPCi(A-(A:), [/■ (fc), <Ji{k))-, 

Submit u\(k) to ISO; 

end 

n 

U\k) = E u\{k) 


end 
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— 'XfXj— Fuses 

3+N+PE 5 conductor cable (3 phases, 
neutral, protective earth) 


Figure 2. Low-Voltage microgrid network structure [15] . 


3.1. System description. The system considered in this case study is composed of 8 buildings 
6 single family houses (SFHs) and 2 multi-family houses (MFHs) connected to a small network. 
The different DER types located at each grid bus/buildings are presented in Table [Tj The 
simulations are performed on an IEEE benchmark LV network (Figure [2|) proposed by the 
CIGRE research group m- The different building and DER models proposed in mUlIIKlH] 
have been implemented in the following simulations. 
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4. Results and Discussions 

This section presents the results generated by applying the proposed control architecture to 
the case study simulation introduced before. 

4.1. VPP operation — Peak shaving. Figure [3] shows the active power flow at the slack 
bus when the cooperation mechanism is implemented ([3^) and not implemented ([SJa) during 
winter time, considering a standard day-night electricity tariff profile (Table [3]). It is shown 
that the DMPC algorithm is successfully maintaining the aggregated load curve in between the 
predicted bounds with the exception of a few overshoots. These violations are correlated to 
strong prediction errors in the uncontrollable electricity consumption which highly influence the 
performance of MFC. Since a cost based performance index has been considered (Eq. 1), the 
different MPC regulators tend to maximize their power consumption during low electricity tariff 
periods (i.e. night time), hence creating virtual consumption peaks. The totally independent 
control MPC formulation performed without any coordination particularly reflects this undesired 
effect with peaks reaching 27 kWe at the slack bus. In order to face this issue, a global power 
flow constraint of ±15 kWe has been added to the optimization problems to attenuate the load 
spikes (Figure [51 d). Figure [31: finally shows the aggregated state of charge (SoC) of the thermal 
and electrical storage units monitored by each MPC. During low tariff periods, the controllers 
heavily charge the domestic hot water (DHW) tanks and batteries to satisfy local (i.e. comfort 
and operating costs) and global (i.e. load constraints) objectives. Moreover, since the feed-in 
tariff of electricity is lower than the night market price (Tabled, the MPC regulator tries to 
maximize self-consumption and the excess PV power generation is recovered by the electrical 
and thermal storage systems. 

In order to analyze the DMPC performance with respect to seasonal variations. Figure 0] 
presents the active power profile when performing DMPC (01)) and when considering only local 
objectives (0^), during summer time. Regarding the large power generation resulting from the 
different RESs installed, the electrical storage capacity (Figure jH;) is not sufficient to satisfy 

Table 1. Parameters for the distributed energy system (DER). 


Building type 

Heating 

Unit Type 

RES 

Unit 

Size 

SFH 

Heat Pump 

Air-Water 

PV panel 

50 [m^] 

MFH 

Heat Pump 

Air-Water 

- 


4 X SFH 

Boiler 

Natural Gas 

PV panel 

100 [m^] 

SFH 

Heat Pump 

Air-Water 

PV panel 

35 [m^] 

MFH 

Cogeneration 

Natural Gas 

- 



Thermal Storage 

Electrical Storage 


Unit 

Size 

Unit 

Size 

SFH 

Hot Water Tank 

34.8 [kWh] 

Battery Stack 

3 [kWh] 

MFH 

Hot Water Tank 

69.7 [kWh] 

- 


4 X SFH 

- 


Battery Stack 

12 [kWh] 

MFH 

Hot Water Tank 

46.4 [kWh] 

Battery Stack 

2 [kWh] 

SFH 

Hot Water Tank 

69.7 [kWh] 

- 
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Figure 3. MFC problem formulation comparison for 3 typical winter days. 
Active power profile at the grid slack bus without (a) and with (b) MFC coor¬ 
dination, as well as the state of charge for aggregated storage unit (c). 


the peak shaving constraint introduced previously. The power consumption/generation curve 
indeed exceeds the ±15 kWe threshold since thermal storage through spacing heating and the 
hot water tank (that saturates over time) is not available during this period of the year, the 
excess of power must be reinjected into the distribution grid. To address this issue, long-term, 
seasonal storage may be necessary through a larger, centralized system located next to the slack 
bus. 
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Figure 4. MFC problem formulation comparison for 3 typical summer days. 

Active power profile at the grid slack bus without (a) and with (b) MFC coor¬ 
dination, as well as the state of charge for aggregated storage unit (c). 

4.2. 24h-ahead Tariff for Load Shifting. As stated previously, regarding the specific per¬ 
formance index considered (Eq. 1), the different MFC controllers maximize power consumption 
during low tariff periods. Hence, in order to steer electricity imports towards specific periods of 
the day, the DSO might provide the ISO with specific day-ahead market prices prior the com¬ 
putation of the 24 hours load prediction. Since the MFC targets are entirely costs based, the 
local DSO may thus be able to shape the aggregated consumption profile regarding its current 
interest. Figured shows the consumption strategy computed by the distributed control scheme. 
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Table 2. Variable electricity tariff profiles. 


Tariff 


Import price Export price 

[CHF/kWhe] [CHF/kWhe] 


Day-Night 


24 h-ahead 


7 am - 10 pm 0.24 
else 0.13 

5 pm - 11 am 0.21 
else 0.16 


0.1 


0.1 


Violation costs [CHF/kWe] 
Predicted profile Global power 

0.5 0.25 


during spring time, while setting lower import costs during mid-day periods, when power gener¬ 
ation from RESs is peaking (Tabled]). Indeed, the maximum of the aggregated power flow curve 
is shifted towards the desired time frame, thus, providing the DSO with an additional ancillary 
service, that is load shifting. 

Table [3] presents the performance values of both pricing schemes considered in this study 
(Figure |S]). As expected, the mean power flow (at the slack bus) during low tariff periods is 
higher than the one during low price periods. As stated previously, the total bound violations 
are related to prediction errors, mainly for the outer conditions. Indeed, the overshoots occurring 
in the morning time are caused by the change in the minimum comfort temperature within the 
buildings and thus, are linked to the outside temperature. However, the peaks observed in Figure 
fig:ahead a are considerably limited since the controllers used the storage systems which have 
been fully charged during the preceding low tariff period (i.e night). A similar behavior can be 
noticed for the strong undershoot in Figurel^i; regarding previous low price period (i.e mid-day), 
the regulators could not provide sufficient storage capacity in order to recover the unpredicted, 
generated power excess. 

Nevertheless, although no grid operation constraints have been included in the local per¬ 
formance, an a posteriori power flow analysis shows that the voltage and phase angle remain 

Table 3. DMPC performance for day-night and 24h-ahead tariffs during 
spring time. 


Tariff 

Mean power flow [kWe] 

Total 


Low tariff 

High tariff 

[kWhej 

Day-night 

10.76 

1.78 

202 

24 h-ahead 

7.88 

2.77 

181.3 


Bound violation [kWhej 

Relative [-] 


Total 

Mean 

Total 

Day-night 

15.6 

1.43 

0.077 

24 h-ahead 

34.4 

4.82 

0.19 
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Hour of year [hr] 


Power profile at slack Bus 
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Figure 5. MFC problem formulation with a modified tariff profile (Tabled]). 

Active power profile at the grid slack bus with the standard day-night (a) and 
the 24h-ahead (b) tariff. 

within a standard quality limit of ±5 p.u. (Tabled]). In order to account for the fast dynamics 
of electrical DERs and hence, their optimal operation, a second control layer should indeed be 
considered during future development. As proposed in [14] for systems with different dynamics, 
a hierarchical approach can be added to the actual MFC problem formulation. The second con¬ 
trol level evaluates the set points of electrical energy systems with a smaller sample time (5-15 
minutes) in order to maintain the grid quality within the desired bounds (±5 p.u.). However, 
to cope with dynamics at higher frequencies (e.g. solar irradiance variations), real-time control 
must still be implemented, which relies outside of the scope of this paper. Curious reader is 
referred to [TO] . 


Table 4. Maximal voltage and angle deviation during DMFC operation. 


Period 

Voltage dev. 

Angle dev. 


[p.u.] 

[deg] 

3-5 January 

0.01 

2.2 

7-9 July 

0.007 

4.6 
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5. Conclusions 

This paper has presented flexible and modular control architecture to define the optimal 
operation plans of DERs connected to a microgrid. Based on an iterative, independent and fully 
connected MFC problem formulation, the presented approach enables VPP operation through 
a central ISO, while providing ancillary services to the local DSO. The case study simulation 
results showed the performance of the proposed control scheme; although a few overshoots have 
been noticed particularly due to strong prediction errors, the coordinating MFC regulators have 
successfully steered towards the desired operation plan. In addition, the latter simulation showed 
the flexibility of the DMPC to steer the aggregated power profile towards specific periods through 
day-ahead market tariffs defined by the DSO. 

Future studies are planned to investigate the optimal design of the energy systems installed at 
each building with respect to the presented case study which has not been shown in this paper. 
Moreover, pricing strategies related to transgressions of the forecasted load curve might also be 
proposed. In fact, in case of constraint violation, the resulting penalty is preserved by the ISO, 
however the latter entity solely represents the aggregated interest of the different DERs and the 
DSO without providing any services or operating strategies. 
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